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ANALYSIS  OF  A  STRUCTURE  WITH  A  RANDOM  GEOMETRY* 

by 

J.  Lyell  Sandera  Jr. 


Introduction 

A  natal  polycryatal,  a  nonwoven  fabric  auch  aa  fait,  a  auapenalon  of 
particlaa  of  randon  aiza,  ahape,  and  location  in  a  vlacoua  fluid,  ara  example  a 
of  what  we  Man  by  a  atructura  with  a  randon  gaoMtry.  Tha  problan  of  intaraat 
in  connaction  with  auch  atructuraa  is  tha  determination  of  aoaa  average  property 
auch  aa  the  effective  Young' a  nodulua  or  the  effective  viacoaity  or  perhepa 
the  probability  denaity  function  of  theae  propart lea.  Solutiona  aetiafactory 
for  practical  purpoaea  to  a  fair  number  of  auch  problena  have  eppeered  in  the 
lltereture,  but  the  author  knowa  of  no  inatance  in  which  aona  (nontrivial) 
average  property  of  a  random  atructura  haa  been  determined  exactly.  The 
difficulty  la  not  in  aolving  aome  known  MtheMtlcal  equation,  but  rather  in 
reducing  auch  problena  to  exact  MtheMtlcal  terma.  The  purpoaa  of  the  praaent 
paper  la  to  becoM  better  acquainted  with  the  nature  of  the  problem  by  Mena 
of  an  exact  analyaie  of  a  very  alnple  random  atructura. 

The  paper  la  concerned  with  the  analyala  of  a  atructura  which  la  eaaentlally 
one-dlMnalonal.  It  turna  out  that  the  analyala  of  thla  atructura  la  much  like 
the  analyala  of  a  Markov  proceaa.  The  fact  that  a  apaclal  coordinate  replecea 
tha  time-like  Independent  variable  of  a  Markov  proceaa  Mkea  for  aome 
complication  aa  might  be  expected  aince  boundary  value  problena  are  generally 
more  complicated  than  initial  value  problena.  The  raaulta  obtained  give  no 
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clear  indication  of  how  to  approach  tha  exact  analysis  of  any  phyalcal  problem 
bat  suggest  that  tha  next  atap  should  probably  be  to  invent  soma  two-dimensional 
problem  sufficiently  simple  to  be  analysed  and  not  so  complicated  as  a  real 
physical  problem  la  likaly  to  be. 


Analysis  of  a  Random  Structure 

Tha  structure  to  be  considarad  is  composed  of  a  uniform  elastic  rod 
pin-jointed  to  a  series  of  elastic  beams  which  are  built-in  to  a  rigid 
foundation  as  shown  in  Figure  1. 
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A— rigid  foundation 
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Figura  1 

The  beams  are  randomly  apaced  with,  on  tha  average,  S  beams  per  unit  length. 

In  other  words  the  probability  of  finding  a  beam  in  a  length  fix  is  equal  to 

Nfix  .  If  the  random  variable  y  la  the  distance  between  successive  beams  than, 
o 

as  is  wall  known  in  such  cases,  the  probability  danslty  function  (p.d.f.)  of  y 

0  .My 

is  given  by  the  Poisson  law  Ne  7  .  eFor  the  present  suppose  the  rod  is  semi- 

infinite  with  one  end  at  tha  origin.  Let  a  force  F  be  applied  to  the  rod  at 

the  origin  and  let  the  resulting  displacement  of  the  end  be  u  .  Let  a  ■  F/u 

be  called  the  stiffness  of  the  structure.  Consider  the  problem  of  finding  the 
•  o 

p.d.f.  ?(a)  of  the  random  variable  a  .  First  add  on  a  beam  at  the  origin. 

The  stiffness  of  the  aupanted  structure  o'  la 
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where  a  is  tha  stiffnass  of  a  beta  with  raspact  to  a  load  at  its  tip.  Tha 
p.d.f.  a'  is  thus 


Q(ef')  -  P(ff)  gr  -  P(o’-a) 

Now  add  a  piaca  of  rod  of  length  y  .  Tha  stiffnass  0"  of  tha  augmented 
structura  is  now  givan  by 

i-  +  _L 

0"  0'  +  K/y 

whara  K/L  is  tha  stiffness  of  a  rod  of  length  L  .  The  conditional 
probability  density  function  (c. p.d.f.)  of  0"  given  0'  is 

T(o"  J0 ' )  ■  R(y)  -  R(Jr  -*7)  — 

whara 

R(y)  •  Na'Ny 

-  0 

and  tha  p.d.f.  of  0"  is 

J T(0";0,)Q(0,)d01 


(2) 


(3) 


(4) 


(5) 


(6) 


Now  tha  augmented  structura  with  the  baasi  and  rod  addad  is  statistically  tha 
same  as  tha  original  structura,  therefore  tha  p.d.f.  of  0"  aust  be  tha  same 
as  tha  p.d.f.  of  0  and  hanca 


P(ff)  ■  ^  / 


P<0'-a)d0' 


(7) 


The  integral  aquation  (7)  la  easily  reducible  to  the  following  differential- 

'O' 

diffaranca  aquation 


(0  P) '  -  XNP  -  -HNP(0-a) 


(8) 


O' 


where  e  prime  denote*  differentiation  with  reapect  to  a  .  If  the  beams  have 
a  random  stiffness  p  with  p.d.f.  S(p)  we  obtain  instead  of  (8)  the  following 


<02P)'  -  KNP  -  -XN  J  P(?)8(ff-?)d? 


(9) 


Equations  (8)  or  (9)  must  be  solved  subject  to  the  normalisation  condition 


J  P(or)dO  -  1 


(10) 


Let  PL(p)  denote  the  Laplace  transform  of  ?(e)  .  We  have 

m  m 

PL(P)  ■  J  P (0)e"epdff  P£(0)  -  -  J  ffP(0)de  - 


PL(0)  -  1 


PL<°>  "  /  P2p(0)<»0  •  V 


(ID 


where  a  bar  over  a  quantity  indicates  its  mean  value.  Taking  the  Laplace 
transform  of  Iq.  (9)  leads  to 

1_SL 

P"  -  KN - -  P  -0  (12) 

L  P  L 

In  the  special  case  8(p)  -  ^  e  for  which  p  -  a  the  exact  solution  for 

P(0)  i* 

P(5> .  Vgm.  ,‘4* + £ 

where  is  the  modified  Bessel  function  of  the  second  kind. 

It  is  interesting  to  make  some  comparisons  with  the  stiffness  of  a  similar 
structure  with  uniformly  spaced  beams.  Let  the  beams  be  spaced  uniformly  at  a 


(13) 
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distance  1/N  and  hava  a  unifora  stlffnaas  a  ,  than  tha  stiffnaaa  of  tha 
atructura  (with  no  baaa  at  tha  flnita  and)  ia 

2 

9  •  (a  KN  +  J-)1/2  -  J  (14) 

Mow  lat  a  -  0  ,  N  "  and  require  aN  to  remain  flnita,  than 

92  ■*  a  K8  (15) 

From  Bq.  (12)  and  with  0  ■  -S^(0)  -  a 

«2  ■  pi*  ■ « »»  <“> 

*  2  ~~2 

tha  same  raault  aa  in  (15),  but  in  ganaral  ( a )  4  o 


Structura  of  Flnita  Lanath 


Tha  aolution  of  tha  problem  in  the  preceding  aaction  for  a  aaml-infinite 
atructure  dapendad  on  a  apaclal  device  which  will  not  work  for  a  atructure  of 
flnita  length.  Conaidar  now  a  somewhat  more  complicated  problem  for  the 


structure  shown  in  Figure  2. 


Figure  2 

Lat  it  ba  required  to  find  tha  Joint  p.d.f.  R(ff,u,x)  .where  0  and  u  are 
tha  stiffness  and  displacement  at  tha  right  and.  Iha  left  and  of  the  structura 
la  subjected  to  a  load  F  as  shown  and  u  is  taken  positive  to  tha  left.  This 
time  augment  tha  structura  by  an  infinitesimal  place  in  such  a  way  as  to  leave 


the  atatietlcel  character  of  the  structure  undisturbed.  Let  R1(9>u>x+Ax)  be 
the  new  R  on  condition  there  is  no  beam  In  the  added  piece  of  length  6x  , 
and  R^  be  the  new  R  on  condition  there  is  a  be aw  In  the  added  length.  The 
unconditional  new  R  is  given  by 

R(9,u,xf$x)  ■  R1(o,u,x+6x)(l-N6x)  +  R2(s ,u,x+6x)N6x  (17) 


since  wa  do  not  or  do  include  a  bean  in  the  addad  place  with  probabilities 
1  -  Nflx  and  N6x  respectively.  Expressions  for  and  R^  era  found  as 
follows . 

1.  No  beams  addad.  Let  91  and  u1  ba  tha  naw  displacement  and  stiffness. 

These  are  related  to  a  and  u  as  follows 

u'  *  u  u  “ 

or 

1  1  fix 

^ma+r 

According  to  the  calculus  of  probabilities 

R1(9' ,u' ,x+4x)do,du'  ■  R(9,u,x)dedu  (19) 


.Ml... 

k-o'&x 


(18) 


where 


dodu  -  J(S2^7)«te'du’  -  (j^TjJ^to'du*  (20) 

and  where  J  is  the  Jacobian  of  the  change  of  variable.  From  this  it  follows 

r1(9,,u,,*+6x)  -  (21> 

2.  Beam  added .  In  this  csss  the  change  in  stiffness  due  to  tha  addad  rod  may 
bs  neglected  and  we  have 

9'  -  9  +  a  (22) 

The  added  beam  exerts  a  force  a  <i'  on  the  old structure  and  causes  a 
displacement  to  the  right  by  an  amount 
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but 


therefor 


In  this  cate 


u«  .  a-«! 

O 


u'  ■  u  -  u" 


u'  ■  u 


or 


u 


g'u1 

o' -a 


3 


o' -a 


and 

^(cr ',u',x+6x)  -  R(g,u,x)  -77^  -  R(er'-a,  ,x)  ^7^ 

Print*  can  now  be  dropped  with  no  loan  in  eenae  and  (17)  becomes 

R(<r,u,x+4x)  -  ,u,x)(“j^)2(l-N6x>  +  ^  R(«r-a,  ,x)N6x 

Subtract  R(a,u,x)  from  both  aide*,  divide  by  6x  and  let  bx  -  0  ,  the 
result  is 


(23) 

(24) 

(24) 

(25) 

(26) 

(27) 


i  I;  to2»)  +  m 


Ha_ 

a-a 


R  (a-a, 


OS_ 

«y-a 


>*) 


If  the  beam  stiffness  is  random  with  p.d.f .  S(p)  than  the  result  is 


sis;!?  (•*»)+« 


9 

»  J  R(<t-?>  ,x)8(5)  ^ 

0 


(28) 


(29) 


These  equations  are  to  be  solved  subject  to  the  usual  normalising  condition 
and  the  initial  condition  R(o,u,0)  ■  0  for  8u  j  F  .  The  equation  for  the 
marginal  distribution  P(c,x)  may  be  obtained  from  (29)  by  integration  over 
all  values  of  u  .  The  result  is 


If  "  X  I?  (®2p)  "  W  +  N  /  p(ff-?>8(?)<i? 

0 


(30) 
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Appropriate  conditions  on  th«  solution  of  this  equation  srs  P(ff ,0)  glvsn  and 

OD 

j  P (o,x)dff  "  1  .  If  ths  left  end  Is  attached  to  a  spring  of  stiffness  or*  , 

0 

for  example,  the  initial  condition  is  P(s,0)  ■  6 (0*9*)  .  Ho  solutions  of 
Eq.  (30)  have  bean  obtained  by  the  author  except  the  one  given  In  the  last 
section  for  the  seal-infinite  structure  when  Iq.  (30)  reduces  to  Eq.  (9)  because 


Properties  at  an  Interior  Point 

Let  ue  now  turn  our  attention  to  the  problem  of  determining  the  joint  p.d.f. 
of  a  ,  u  ,  and  force  f  (In  the  rod)  at  a  point  x  Interior  to  a  structure  of 
finite  length  L  .  This  problem  may  be  solved  In  terms  of  previously  determined 


0 


Figure  3  x  L 


probability  density  functions  if  we  make  a  cut  In  the  structure  at  the  point  x 


and  then  determine  the  forcas  and  displacements  necessary  to  restore  continuity 


Figure  4 
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Flgure  4  >how*  the  cut  structure  near  the  point  x  .  Let: 
a,  ■  the  etlffneae  et  the  point  1 

9^  “  the  stlffnese  et  the  point  2 

a  ■  the  etiffneea  of  the  continuous  structure  et  x 

f  -  the  force  necessary  to  effect  the  closure  which  equals  the  force  et 

x  on  the  continuous  structure 

Uq  “  the  gap  between  1  end  2  before  the  two  pieces  of  the  structure 
ere  joined  together 

u  -  the  dlsplecesMnt  et  x  of  the  continuous  structure  (positive  to  the 
left) 


u^  -  displacement  of  point  1  after  the  force  f  Is  applied 
u2  -  displacement  of  point  2  after  the  force  f  Is  applied 

We  have 


«1  “  uo  •  f/°! 

(31) 

U2  "  f/02 

(32) 

The  structure  will  be  continuous  and  In  equilibria  If  f  and  u  are 

determined  from 

U1  “  u2  -  u 

(33) 

This  gives 

f  fflVo 

tTl+°2 

(34) 

Vo 

al*2 

(35) 

The  stiffness  of 

the  joined  structure  Is  obviously 

a  m 

(36) 
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Now  tinea  (?,  and  uQ  ara  independent  of  we  may  write 

u0,ffl'®2 

p(«T,u,f,x)  ■  R(g x) P (q2,L-x)J ) 
where  p  la  the  required  p.d.f.  and  R  and  P  have  been  prevloualy  defined. 


(37) 


Bquatlona  (34),  (35)  and  (36)  may  be  solved  for  and  uQ  with  the  result 


uQ  -  eru  /(ffu-f) 
8^  m  o  -  f/u 
<?2  -  f/u 

Fro*  (38)  we  calculate  J  to  be 

J  »  8 / (ou-f ) 

The  result  Is  thus 


(38) 


(39) 


p(8,u,f,x)  -  Ha  -  l  ,  >*)?(„ 


8  >  “  >  0 
u 


-  0 


-  <  0  or  8  <  - 
u  u 


(40) 


With  a  sufficient  amount  of  effort  one  can  derive  an  equation  resembling 

a  Fokker-Planck  equation  satisfied  by  p  ,  but  the  results  for  L  finite  or 

Infinite  differ  from  each  other.  The  author  has  been  unable  to  arrive  at  these 

equations  from  first  principles  or  to  determine  whether  or  not  these  equations 

together  with  certain  boundary  conditions  constitute  an  appropriate  formulation 

2 

of  the  problem.  The  equation  In  question  contains  but  not  and  the 

X  dx* 

association  of  such  an  equation  with  two  point  boundary  conditions  seems 
Inappropriate.  On  the  other  hand  there  does  not  eeem  to  be  any  reasonable 
second-order  Fokker-Planck  equation  satisfied  by  p  . 
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